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ABSTRACT 


This thesis addresses the evaluation of the apparent target-background temperature 
difference (ATapp) at maximum range and compares ATapp with MDTD/MRTD of a typical 
FLIR system. The atmospheric propagation code (SEARAD) and Planck’s radiation law 
were employed to obtain atmospheric transmittance and path radiance. The atmospheric 
parameters were selected as close to the aircraft overflights as possible using PREOS 92 
data as inputs of the SEARAD code. Ship target (RA^ POINT SUR) modeling was 
established for MDTD/MRTD evaluation using a rectangular parallelepiped model of the 
ship’s physical length, width and height. The geometry data for MDTD/MRTD 
evaluations were also selected from the PREOS 92 experiment measurement set. The 
MDTD/MRTD functions for a generic FLIR in wide field of view (WFOV) application 
were deduced from Shumaker. Johnson criterion was employed as a detection criterion. 
Resolution line-pairs at detection range to resolve the target have also been evaluated and 
compared against Johnson Criterion. 

The temperature differences between AT^pp and MDTD at detection range show 
large scatter, ranging from 5% to 600 %. They also show agreement with the same sensor 
altitude and viewing angle. A comparison of ATapp with MRTD at classification range 
and identification range show that using NFOV would be more appropriate for target 
classification/recognition. Resolution line-pairs at detection range derived from a typical 
WFOV MRTD curve show 10%~50% error for the eight runs, which is acceptable as 
compared against Johnson Criterion quoted as 1+ 0.25 mRad. 
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I. INTRODUCTION 


Infrared (IR) systems have been under development continually over the past 
decade. IR systems are usually used for target search, detection, classification and 
identification. The performance of many military systems using FLIR technology is 
strongly dependent on environmental conditions in the tactical operation area. The crucial 
performance parameter of a FLIR system sensor is its maximum detection range (MDR) 
for a specific target. MDR is determined by application of a defined apparent temperature 
difference {AT) criterion for detection or recognition/classification at a selected 
probability. A generic tactical scenario situation is shown in Figure 1 when an airborne 
FLIR sensor detects a ship at a certain distance. In order to predict the MDR, it is 
necessary to determine the target minus the background temperature contrast {AT) as 
viewed by the sensor. To be detected the point target AT must equal the Minimum 
Detectable Temperature Difference (MDTD) or Minimum Resolvable Temperature 
Difference (MRTD) of the FLIR system. The temperature contrast calculation requires an 
accurate prediction of the spectral radiance of the target and the background, and the 
attenuation and radiance of the path between the sensor and the target. From experience 
we know that the meteorological conditions have a great influence on the performance of 
FLIR systems. Over the past several decades several atmospheric propagation codes 
(HITRAN, MODTRAN and LOWTRAN) were developed by the Phillips Laboratory 
for both high resolution and low resolution transmission. The HITRAN codes, for 
example, FASCODE , have very high resolution and resolve individual molecular rotation 
, vibration lines. But they are slow and expensive. LOWTRAN (the latest version is 
LOWTRAN 7) is essentially a curve fitting match to empirical data. It was designed to 
calculate atmospheric transmittance and radiance averaged over a 20 inverse centimeter 
frequency interval with steps at multiples of 5 inverse centimeters, and accuracy about 10 
%. This is a faster code and has received wide acceptance. MODTRAN is an atmospheric 
propagation code that can be used at moderate resolution, to 2 cm'*, using data stored at 1 
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cm intervals. M0DTRAN2 also extends the spectral range to include the near ultraviolet 
through the microwave. Recently, a new research atmospheric propagation code called 
."SEARAD" (modified MODTRAN2) has also been developed to predict the 
transmittance, path radiance, and sea and reflected sky radiances [Ref 1]. SEARAD 
code adds sea surface radiance based on Cox-Munk (Cox and Munk, 1954, 1956) 
statistical model [Ref 1] which also includes sun glint effects. It is a DOS-compatible 
program which can run on a personal computer and calculates sea radiance. 

This thesis will address the calculation of the apparent AT zA maximum range and 
compare AT with MDTD/MRTD to examine the influence of the meteorological 
parameters that cause a variance in FLIR performance. The source of meteorological 
parameters was based on the PREOS92 experiment that was conducted on August 04, 
1992. The PREOS computer program [Ref 2] is a NRaD developmental Tactical 
Decision Aid (TDA) model that was intended as a range prediction technique for airborne 
FLIR systems. PREOS92 represents an experiment for validation of the PREOS program 
conducted off the California coast from Monterey in 1992. Details of this experiment will 
be described in Chapter III. In this experiment. The research vessel POINT SUR (R/V 
POINT SUR) was used as a target and as a measurement platform. The skin temperature 
distribution was recorded throughout the experiment using a set of 15 thermistors on the 
ship skin. Thermal image data obtained from the R/V POINT SUR measured from navy 
aircraft were used with the ship temperature data to determine the average ship 
temperature. Detection range, classification range and identification range were 
determined by military aircraft; overflights. Meteorological and sea surface data measured 
on board the POINT SUR were used as inputs to the atmospheric propagation 
computations. A Radon counting instrument on the R/V POINT SUR was used to 
determine the air mass parameter (ICSTL) representing the degree of coastal aerosol 
influence for inclusion in the SEARAD propagation code. In this experiment, the NRaD 
aircraft, the military aircraft and the radiosonde balloons could not occupy the air space 
around the R/V POINT SUR at the same time. The experiment was planned to have data 
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from those sources as closely spaced in time as possible. The meteorological parameters, 
observed detection ranges and aircraft altitudes selected in this experiment were used as 
inputs to the "SEARAD" code to calculate the atmospheric infrared (8 to 12 p,m) 
transmittance and path radiance. 

The objective of the calculation of apparent temperature difference {AT) was to 
evaluate the effective Minimum Detectable Temperature Difference (MDTD) and 
Minimum Resolvable Temperature Difference (MRTD) for the FLIR system from the 
PREOS 92 measurement set. The Johnson criterion, the number of line pairs need across a 
target, was used for evaluating the MDTD/MRTD of the FLIR system. The /17 calculation 
was based on the meteorological data collected that was closest in time to the aircraft 
overflights on the test. After accessing the transmittance and path radiance calculated from 
SEARAD, Planck's radiation law was applied to find the apparent radiance difference and 
corresponding temperature difference {AT) at the detection range. Generic FLIR 
parameters were used to estimate MDTD/MRTD of the FLIR sensor. Furthermore, the 
number of line pairs required to satisfy the MRTD criterion at the detection range were 
investigated to improve the performance of prediction. 
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11. THEORY AND BACKGROUND 


A. DETECTION SCENARIO 

Figure 1 shows a typical scenario to detect a ship with an airborne FLIR. [Ref. 2 ] 
In order to model the detection of the ship, we have to know or calculate the detected 
radiance seen by the sensor, after propagation through the atmosphere. From this we can 
compute an apparent temperature difference AT from the ship and background and 
compare this to the Minimum Detectable Temperature Difference (MDTD) or Minimum 
Resolvable Temperature Difference (MRTD). If the apparent temperature difference is 
equal to the MDTD or MRTD, the target should be detected or resolved. 

In order to find the apparent temperature differences between ship and 
background, the ship and the background radiance must be adjusted for the atmospheric 
effects over the spectrum of the sensor, and then the effective temperature difference 
calculated. The ship radiance must be adjusted by the transmittance and path radiance 
effects. Also, the background flux has to be adjusted for path radiance, sea surface 
radiance and reflected sky radiance. [Ref 2, p. 3] 

B. RADIATION LAW 


The following thermal radiation laws were used for the IR radiation and apparent 
temperature difference calculations. These equations are in general use and are 
summarized in Ref 3 and Ref 4. They are reproduced here for convenience. 

1. Planck's Blackbody Radiation Law 

A blackbody is an idealized source of radiant energy. The emitted energy depends 
on the temperature of the object and the wavelength. The spectral radiant emittance 
from a blackbody is given by Planck's blackbody law and is shown in Equation (1), 





f . 1 

W' J 

1 

_ y 


(1) 
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where 


X = wavelength (m), 

k = Boltzmann's constant = 1.38x10'^^ (W - sec/K), 
h = Planck's constant = 6.626x10'^'* (W-sec^), 

T = temperature (K). 

Planck's law can be applied to all FLIR systems and is valid for the entire electromagnetic 
spectrum. The spectral radiant emittance of a blackbody ranging from 500 to 900 K is 
shown in Figure 2. [Ref 3, p. 55] From this figure it shows that the spectral radiant 
emittance dependents on the temperature of the object and the wavelength. 

1 . Stefan-Boltzmann Law 

Integrating M(X,T) over all wavelengths gives the total radiant exitance 
Equation (2), from a blackbody surface in terms of temperature T. 

( 2 ) 

where 

a = Stefan-Boltzmann constant = 5.6697x10-* (W/m^/K'^) 

For a non-blackbody source, this law is modified by the effective source emissivity 8, 
giving the form in Equation (3) 

M =£(tT\ (3) 

For in-band radiance, the radiance N can be calculated using the universal blackbody 
curve as shown in Figure 3 or the proportional radiation table (Appendix A) with Equation 
(4) [Ref 5]. 

saP 

N = - Aq , (4) 

n 

where 

s = emissivity of the source. 

Aq = difference of two q values from the proportional radiation table (Appendix A) 
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The proportional radiation q - f (A,T) represents the fraction of the radiant exitance 
emitted by a blackbody at temperature T at all wavelengths up to the selected value of A. 
It can be described in Equation (5): 

l‘u(X,T)dA 

If q values are not shown on the table, they can be calculated by interpolation between 
two q value entries. 

3. Lambert-Beer Law 

The atmospheric transmittance is a strong function of wavelength, which can be 
described by the Lambert-Beer law. 

ri^) = e-^ , ( 6 ) 

where 

T (A J = spectral atmospheric transmittance, 

A = wavelength, 

H (t) = extinction coefficient, 

R = path length. 

The extinction coefficient depends on the atmospheric molecular composition and 
aerosol concentration. It is very sensitive to the specific frequencies of the molecular 
vibration and rotation transitions. 

4. Absorption and Scattering 

The infrared radiation will attenuate as it propagates through the atmosphere by 
the processes of absorption and scattering. The attenuation is characterized by the total 
extinction coefficient. 

= (7) 

where 


7 



// = total extinction coefficient, 

= extinction coefficient for total absorption, 

//g = extinction coefficient for total scattering. 

Both the scattering and absorption coefficients have two components; one fi'om 
the molecules of the air and the other fi-om aerosol particles suspended in it. Scattering by 
aerosols and molecules has a large effect in the visible region, while absorption dominates 
in the infi-ared region of the electromagnetic spectrum. Since FLIR systems are designed 
to operate between 8 to 12 pm, it is important to predict the transmittance of the 
atmosphere as a function of wavelength and weather conditions. These two coefficients 
are shown in Equation (8) and (9). 

( 8 ) 

Ms=<^m + cr^ (9) 

where 

k^= molecular absorption coefficient, 

A: aerosol absorption coefficient, 

cr^=molecular scattering coefficient, 
cr^=aerosol scattering coefficient. 

The relative values of the four coefficients depend strongly on the density and molecular 
composition of the atmosphere and the composition and size distribution of the aerosols. 
[Ref 6, p. 12-7] 

A range of computer models have been developed for computation of atmospheric 
transmittance for various defined meteorological scenarios. These models include 
LOWTRAN, MODTRAN and FASCODE. Recent variations of these codes, for example 
SEARAD, also compute path radiance, and the background radiance can be converted 
into background temperature directly. [Ref 1] 

C. RADIATION IN TDA MODELING 

A Tactical Decision Aid (TDA) code such as a FLIR range prediction generally 
consists of three models: 
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• A target model which describes a target as an average temperature difference 
A T against a background. 

• An atmosphere model which describes the radiation in apparent temperature 
difference from the target to the sensor. 

• A detection model describes the sensor performance for specifying MDTD 
for detection or MRTD for classification or recognition. [Ref 7] 

The PREOS model predicts the maximum detection range of an airborne FLIR system 
over the 8 to 12 micron band. [Ref 2] In this model, the SHttPSIG model provides the 
ship target signature based on ship skin temperatures and meteorological data. The target 
signature appears as an average temperature between the ship and the background. A 
modified form of the LOWTRAN6 model computes the propagation of the radiation 
difference to the detection range. A parameterized MDTD equation based on a known 
system performance represents the FLIR characteristics. 

The following equations and definitions to describe TDA modeling are taken from 
C.P. McGrath [Ref 2, pp. 3-6], Jon C. Kreitz [Ref 3, p. 12-15] and C.R.Zeisse [Ref 1, 
P-7] 

1. Target Model 

Typically a target model consists of an average ship temperature and the difference 
fi-om the adjacent background temperature AT^. To incorporate atmospheric changes in 
this quantity it must be converted into an equivalent radiance difference. This is done 
through an integration over Planck's blackbody distribution over the system bandwidth 
and gives a zero range equivalent in-band blackbody radiance N(H, $), for the target, 
where H represents the altitude of a FLIR sensor and 6 represents the elevation angle 
from the target to the sensor. An equivalent radiance for the background must be 
included. 
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2. Atmosphere Model 

The atmosphere model computes the path transmittance, the path radiance (e.g 
from SEARAD) and the background radiance N(H, 0 )^. N(H, is the sum of the three 
radiance terms: 

• Sea Surface Radiance N(H, 6)ss ; 

• Reflected Sky Radiance N(H, ; 

• Path Radiance N(H, 0)^. 

These three terms may be calculated from the modified MODTRAN2 (SEARAD) code, 
and can be described by Equation (10). 

n{h.0 \ = N(H,e + N{H,e + n{h.0 ), ( 10 ) 

Also the radiance of the background has been converted directly to an equivalent 
background temperature by the SEARAD code. [Ref 1] 

3. Detection Model 

The detection model compares the apparent temperature difference A 7^, (H, 0 ), 
which varies with the range and viewing angle, to the threshold parameter of the sensor. 
AT^^j^(H, 0 ) is the apparent target minus background temperature difference seen by an 

infrared sensor, and is a function of the sensor's altitude H and viewing angle 0 relative to 
the zenith. It can be described by Equation (11) 

where 

t[h,0)^^^ = average apparent target temperature, 

T[h,0 = average background temperature. 

The major purpose of this calculation is to deduce the apparent temperature difference 
A Tapp at the observed detection range, and compare this with the MDTD/MRTD of the 
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aircraft FLIR sensor, to evaluate the effectiveness of the MDTD/MRTD criteria for the 
FLIR system. 

D. DETECTION CRITERION 

1. Introduction 

The target characteristics play an important role in a FLER system analysis. The 
primary reason for describing the properties of targets is to provide a basic detection 
criterion for analyzing the interaction of the target and the FLIR system. Two parameters 
frequently used to characterize targets for purposes of calculating their interaction with 
FLIR systems are the “critical dimension” and “average effective temperature difference.” 
An “aperiodic” detection criterion is sometimes used, meaning that the irradiance at one 
pixel in the scene exceeds a defined threshold to be identified as a target, with no minimum 
number of resolvable line pairs on the image. These target features provides a way to 
deduce the ability of an observer using a FLIR to obtain information regarding scenes, 
objects and targets from their displayed images. However, other analyses believe that 
some level of resolution above one scan line of the target compared with the background 
is necessary for detection. 

2. Target Transform Function 

Laboratory measurements will describe the performance of the FLIR-aided eye in 
the detection of an isolated rectangular temperature target or resolution of the bars in a 
standard “bar pattern”. Rectangular target detection is a form of “aperiodic detection”, 
characterized by a FLIR parameter called the Minimum Detectable Temperature 
Difference (MDTD). Resolving a bar chart represents “periodic detection” and is 
characterized by a FLIR parameter called the Minimum Resolvable Temperature 
Difference (MRTD). [Ref 8, pp. 2.1-2.4] These two parameters will be described in detail 
in the next section. 

Evaluating the field performance of a FLIR requires the selection of an idealized 
target to characterize a particular visual task. This is a judgmental problem for which no 
firm rules are established. In general, detection of a real target is equivalent to detection 
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of an isolated rectangle whose area is equal to the target area and whose contrast 
temperature difference AT to the background is equal to the average AT between the 
actual target and background. Another visual discrimination capability involves resolving 
bars of a bar pattern whose AT between hot bars and cold bars is set equal to the average 
AT between the target and background, and whose spatial frequency is determined by the 
target's critical dimension and the task selected. It is helpful to represent real targets as 
idealized MDTD and MRTD test targets whose detectability/resolvability can be predicted 
and measured in the laboratory. [Ref 8, p. 2-5] 

3. Equivalent Bar Chart Selection 

Several criteria exist to relate the resolution of bar charts to levels of visual 
discrimination. The first one is the “Johnson Criterion,” which has been widely used for 
land vehicle targets. [Ref 8] In 1970, Johnson first put a set of tribar patterns next to a 
military target under different conditions. He used the minimum dimension of the target as 
a reference. [Ref 8, p. 2-6] But Moser found that the Johnson Criterion was not suitable 
to ship targets. The Moser study involved many observations of cathode ray screen images 
of ship silhouettes formed of black squares of varied scale (representing “pixels”). These 
studies produced the following standard numbers of pixels needed for : detection (1 
pixel), classification (66 pixels) and recognition (400 pixels). Shumaker relates the pixel 
criteria advanced by Moser to the equivalent bar chart spatial frequency. Moser found that 
the minimum ship dimension had little significance, but classification and identification of 
ships was related to the number of pixels falling within the area of a ship image. Table 1 
shows the data relating the Moser and O’Neill conclusions to the Johnson criteria for 
detection, classification and recognition. Moser also found that the square root of 
projected area was “critical dimension” that gave consistent predictions for very elongated 
targets, while the “minimum dimension” led to incorrect conclusions. In lieu of prior target 
information, Moser recommended that this “average dimension” be used. [Ref 8, pp. 2.6- 
2 . 10 ] 
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4. Determination of Target Spatial Frequency 

In FLIR performance analysis, the effective target area is the projected area normal 
to the line of sight from the sensor to the target. The area dependents on the target 
dimensions, the observing aspect and the obscuration of the target. The physical 
dimensions of a target must be known prior to performing radiometric calculations. 
Representing the target as a rectangular parallelepiped of dimension of a target’s length, 
width and height simplifies subsequent calculations. 

Another important consideration for FLIR analysis is the target aspect, the 
elevation angle 6 and azimuth angle (j) seen from the front view. The area of a target 
image as a function of its dimension and the azimuth and elevation angle is 

Aj =lws\n 9+1m cos 9 s\n (j) + hi cos 9 cos <!> , (12) 

where 

h = the unobscured target height, 
w = the target width, 

(p = the target length, 

/ = the azimuth angle, 

9 = the elevation angle. 

Figure 4 shows the orientation of a ship target model for calculating the projected area. 

Using Moser’s suggestions the critical dimension is then =yf^. Dividing 
the critical dimension by the number of bars N in the bar chart gives the physical 
dimension % to be resolved to accomplish the task, i.e. 

Finally, the spatial frequency characterizing the target can be written as lines per 
milliradian 

= (cyc/raRad) (13) 

where is the range in meters from sensor to target. [Ref 8, pp. 10-16] The relation of 
the spatial frequency to range R and physical dimension x are shown in Figure 5. 
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THE FLIR SYSTEM PERFORMANCE 


E. 

The purpose of a FLIR system is to detect the arrival of IR energy, electronically 
■process the detected energy, and display the results to the FLIR operator. The ability of 
the FLIR system to perform its function depends on how well it will provide a scene that 
the operator can view directly. 

Two standard descriptions of FLIR performance are : 

• Minimum Detectable Temperature Difference (MDTD) 

• Minimum Resolvable Temperature Difference (MRTD) 

Each of these provides a sensitivity measurement of the FLIR system. Both of them 
employ the concept of the Equivalent Blackbody Temperature of a non-blackbody object. 
The equivalent blackbody temperature is that temperature a blackbody must have to emit 
the same radiance as the target. This allows the use of a single value to represent the 
target's temperature and emissivity. 

1. MRTD Description 

MRTD is a function rather than a simple number. It gives the AT between the hot 
and ambient temperature bars of a standard 4 bar (7:1 aspect ratio) chart required to make 
the bars just resolvable as a function of the spatial frequency of the bars. MRTD is a 
measure of the performance of the entire FLIR system. Since the MRTD includes the 
observer, the measurement is subjective. [Ref 8, p. 8-52] The MRTD can be written in 
terms of the Noise Equivalent Temperature Difference (NETD) and Modulation Transfer 
Function (MTF) of the FLIR system, the detector dwell time on a target point, and the 
required signal to noise ratio, but must also include a representation of the human operator 
perception factors which may be modeled as an eye filter function and integration time. 
[Ref 8, pp. 8.37-8.38] It typically shows a rapid increase to an asymptotic value at a 
spatial frequency beyond which the pattern can not be resolved at any temperature 
difference. For a typical FLER (Common Module), the instantaneous field of view (IFOV) 
is approximately 0.25 mRad so that the in-scan detector subtense Ax is approximately 0.25 
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1 1 

mRad and the cut off spatial frequency (asymptote) is approximately — =-= 4 

cyc/mRad. If we use the “rule of thumb” asymptote between 0.7/zlx and 0.9/zlx , we 
would find that the asymptote is approximately 3 cyc/mRad. The effective IFOV depends 
on the magnification of the optics which may be variable in some FLIRs. 

Shumaker gives a typical form for the MRTD. [Ref 8, p. 8-52 ] 

^ ^ mtf3(v) L l J 

where SNRT = signal to noise ratio threshold, 

NET = noise equivalent temperature , 

Px = noise filter factor, 

Ax = in-scan detector subtense (mRad), 

Ay = cross-scan detector subtense (mRad), 

MTF = modulation transfer function, 

v = spatial frequency for which MRT is desired in cyc/mRad, 

L = length-to- width ratio of the bar ( 7 ), 
te = eye integration time in second, 

Fr = frame rate s'\ 

Nos = overscan ratio, 

Nss = serial scan ratio. 

Shumaker also provides an expression for the Noise Equivalent Temperature 
Difference (NETd), [Ref 8, pp. 8-7] 


NETj,- 


io(fov. fov^ f, n„ n„) 




where FOVx = in-scan FOV in mRad, 

FOVy = cross-scan FOV in mRad, 


= number of detectors. 
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r| sc = scan efficiency, 

5N/5T = derivative of Planck’s equation in W/cm^ /K/Sr, 

D** = the band average detectivity in cm AV, 

D = aperture diameter in meters. 

To = transmission of the optics. 

For calculation purposes it is convenient to write MRTD as function of spatial 
frequency f in a parametric form. 



r 

f \ 

/ 1 X 

MRTD^SL + 

SC 

2/ 

1 

exp[^;£J?V^J 


\ 


J 


where SL (sensitivity limit), SC (sensitivity constant) and ER (Equivalent resolution) are 
FLIR constants derived for any given FLIR by fitting the equation above to the measured 
or calculated MRTD of the system for a 4 bar target pattern with ratio ^=7 of bar length 
to width. [Ref 9] FLIR constants are not generally available for specific FLIR instruments 
but can be deduced from suitable generic models. For example, Shumaker example 8-13 
shows the typical FLIR parameters to model a generic FLIR system. [Ref 8, p. 8-59] 
Figure 10 shows the relation between the MRTD and spatial frequency for a typical FLIR 
system. Since MRTD is a function of spatial frequency, the maximum range for 
classification or recognition is the range when the ship/background AT matches MRTD at 
the spatial frequency given by the appropriate criterion, (e.g., the Johnson criterion) 

2. MDTD Description 

The MDTD definition of a FLIR system is similar to the MRTD , but the target is 
represented as a large square and MDTD is not as dependent on the eye parameters. The 
definition of MDTD is the minimum detectable temperature difference between a large 
square target and its adjacent background that just makes the square target detectable. 
[Ref 8, pp. 8.66-8.67] The square target represents a target of known angular subtense 
and temperature difference to the FLIR system. Also, the MDTD of a FLIR system can be 
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tabulated as a function of target angular subtense (TAS) for table look-up purposes. [Ref 
10, p.13-15] 

Shumaker provides a typical form for MDTD 


^J.(v) = (SNRT)(NET)(n,+r-)(A.A,.) 

'{rl+rl+Slr)‘.F,N„N. 


4 


(17) 


where Qt = target angular subtense (mRad^), 
rs = system resolution (mRad), 
re = back-end resolution (mRad). 

The other definitions of the parameters in Equation (17) such as SNRT, NET, Ax, 
Ay, te, Fr, Nos, Nss, are the same as those in Equation (14). 

The calculation of MDTD for a generic FLIR system can be described by a curve 
fit equation as shown in Equation (18). [Ref 9] 

r 1 ^ 


MDTD = SL + 


sc*(tas^ 


+ ER^Y /TAS^ 


(18) 


where SL , SC and ER are empirical FLIR constants as before. TAS represents the target 
angular subtense in mRad which equals 10® Aj/R^. Figure 12 shows the relation between 
MDTD and target angular subtense (TAS) for a typical airborne FLIR system derived 
from equation (17) and equation (18) for WFOV application. This can be converted into a 
function of range for a specific target. 
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III. OVERVIEW OF PREOS 92 EXPERIMENT 


A. INTRODUCTION 

The data for this thesis were obtained from an open ocean field experiment 
conducted from 27 July to 4 August 1992 off the coast of Monterey, California by Naval 
Command, Control and Ocean Surveillance Center (NCCOSC), Research, Development, 
Test and Evaluation Division (NRaD) in cooperation with the Naval Postgraduate School 
(NPS). The purpose of the experiment was to collect a comprehensive set of 
meteorological and ship temperature data, infrared images, and airborne FLIR detection 
range data to evaluate the Performance and Range of Electroptical System (PREOS) 
detection range algorithm developed at NCCOSC. Data and preliminary results of the 
PREOS comparison were documented by Kreitz [Ref 3] and more completely by Hughes 
and McGrath. [Ref 11] 

The Research Vessel POINT SUR (R/V POINT SUR) was the primary target 
platform. R/V POINT SUR is a 135 foot ship (dimension ; 41.5Lx9.75Wx8.8H (m)) 
owned by the National Science Foundation and operated by Moss Landing Marine 
Laboratory. It served as a platform for meteorological and surface measurements and as 
an infrared target for the airborne FLIR and imaging systems. Scientists from NPS, NRaD, 
and Naval Surface Warfare Center (NSWC) were aboard the POINT SUR to conduct the 
onboard measurements. Air temperatures ranged from 11.5 to 17 8° C during the 
measurement. Sea temperatures were from 11.0 to 18.1 <>0. Wind speeds ranged from 3 to 
26 knots. Stratus clouds dominated until the last day of the test when clear skies prevailed. 
The stratus base ranged from 140 m to 300 m according to the aircrew reports. Vertical 
profiles of atmospheric pressures, relative humidities, temperatures, and dew point were 
collected from radiosonde launches from the R/V POINT SUR. Due to the prevailing 
weather conditions the analysis of this thesis is applied only to the data of 4 August. 
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B. DATA MEASUREMENTS 


1. Shipboard Measurements 

Skin temperatures of the ship were collected from 15 thermistors attached to the 
ship, and a hand-held radiometer. The thermistors were mounted on the large area 
surfaces of the ship and were continuously recorded. Spot measurements of infrared skin 
temperature were made several times daily with a hand-held radiometer that operated in 
the 8 to 12 micron band. Meteorological data were collected by the NPS from a variety of 
ship board sensors. The parameters measured included wind speed, wind direction, 
relative humidity and pressure. Atmospheric radon concentrations were also recorded to 
determine the influence of continental versus marine air masses. The ship navigation data 
include latitude, longitude, gyro heading, and ship speed. 

2. Airborne Measurements 

The NRaD airborne platform was a Piper Navajo aircraft equipped with a thermal 
imaging recording system. The aircraft flew 11 missions. A total of 2603 image frames of 
POINT SUR were collected with the AGA 780 Radiometric Thermal Imaging System, 
which was operated in the 8 to 12 micron band. 

3. FLIR Measurements 

The R/V POINT SUR not only served as the primary meteorological measurement 
platform, but also was used as the target for the FLIR and AGA infrared imaging systems. 
Seven naval aircraft each employing an 8 to 12 pm FLIR system flew range detection 
missions during the experiment. Eight missions were flown, four at night, and four during 
the mid-morning, making a total of 57 passes at the target. The measured detection 
ranges, classification ranges and identification ranges during the overflights provide the 
primary foundation for evaluating the temperature difference transmitted to the FLIR 
sensor. During the detection runs, all detection ranges were determined using the wide 
field of view lens for common module FLIR systems with dual FOV. A target detection 
occurs when the aircraft operator first detected a dot on the display. This is a subjective 
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value. Display images were preserved on video tape which is now archived at NPS. [Ref. 

■ 3, 12] 
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IV. DATA SELECTION 


The data utilized for the calculation of apparent temperature difference AT \n this 
thesis were selected from the data gathered in the FLIR92 experiments. The data of the 
last day, 4 August 1992, were selected for use for the reasons outlined by Kreitz [Ref 3], 
namely that the overcast conditions obtaining during the other overflights seriously 
impacted the application of the atmospheric propagation codes. Table 3 shows eight 
aircraft overflights conducted on 4 August, during the time period from 0411 to 0522 
(GMT). Data from a radiosonde launch at 0401 (GMT) provides input information to the 
SEARAD code and are listed in Table 4. Table 5 shows the meteorological data needed 
for atmospheric propagation code (SEARAD). Atmospheric transmission and path 
radiance values obtained from the SEARAD code for calculating the apparent temperature 
difference at detection range, classification range and identification range are listed in 
Table 6, Table 7 and Table 8 separately. Table 3 shows the observed detection, 
classification and recognition ranges for eight flights that were used to evaluate the FLIR 
MDTD and MRTD. The average radiometric ship target temperatures corrected to zero 
range have been analyzed by Hughes and McGrath [Ref 11] and are listed in Table 9. The 
meteorological data were selected to be as close as possible to the aircraft overflight time 
and provide a basis for apparent zlT, MRTD and MDTD calculations. These data in Table 
3, Table 4 and Table 5 served as inputs to the SEARAD code to obtain the transmittance 
and path radiance between the detection range, classification range and identification 
range from the aircraft to R/V POINT SUR. 
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V. ESTIMATION OF SYSTEM MDTD AND MRTD 


To complete the comparison of apparent AT with the system requires numerical 
values of these quantities as functions of spatial frequency for the specific sensors used. 
Since this detailed information could not be obtained for the purposes of this thesis, 
suitable model functions were constructed using the available information. This section 
describes the generation of the comparison performance parameters. 

From the widely distributed manuals on the Marine Patrol Aircraft sensor suites the 
FLIR was expected to be a Common Module type with 180 Mercury Cadmium Teluride 
detectors in parallel scan, with wide field of view (WFOV) 15** elevation and 20® azimuth, 
and 5® elevation and 6 ®40’ azimuth in Narrow Field of View (NFOV). Moser [Ref 13] 
has used detector angular subtense Ax of 0.25 mrad and Noise Equivalent Temperature 
0.3 ®C for such a system. 

FLIRs of the Common Module generation have been discussed extensively by 
Shumaker [Ref 8], Lloyd [Ref 14] and Howe. [Ref 15] 

A. GENERIC FLIR PARAMETERS FOR NFOV 


The theoretical MRTD and MDTD can be expressed in Equation (14) and (17) 
discussed in Chapter II, Section E . [Ref 8, p. 8-55] 


MTF,(v) 


1/ r 


T 


-,v. 


,^eK^OS^Ss] 


'^2 


UDT(v) = (SNRT)(NET)(n,+.,^ )(AxAj.)' 


n. 


+rg +^T)t^F^N q^N 


(14) 


(17) 


Shumaker provides the following set of parameters that are postulated to be appropriate 
to the Common Module FLIR in NFOV mode. [Ref 9, p. 8-59] 


25 


7 


FOVx 

6.86" 

FOVy 

R 

30 sec’’ 

Nos 

Nss 

1 

Nd 

rise 

0.75 

To 

D 

0.10m 

Ax 

Ay 

0.25 mRad 

D* 

f 

0.20 m 

MTF(v) 

te 

0.1 sec 

SNRT 


5.16" L 

1 Te 

180 u 

0.70 Ns 

0.25 mRad NET 


4X10'“cmHz^^ AVatt 
Tabulated below 
2.5 


0.3 mRad 
0.08 
2 

0.1 °K 


MTF(v) calculation for NFOV 


Spatial Frequency (v) 

Ax= 0.25 mRad 

1 mRad 

(cyc/mRad) 

MTF 

MTF 

0.00 

1.00 

1.00 

0.25 

0.89 

0.72 

0.50 

0.79 

0.40 

0.75 

0.64 

0.14 

1.00 

0.53 

0.00 

1.25 

0.43 

0.00 

1.50 

0.35 

0.00 

1.75 

0.27 

0.00 

2.00 

0.20 

0.00 

2.50 

0.10 

0.00 

3.00 

0.04 

0.00 


Px=(I+(2v» is computed from the equivalent back-end resolution which is 

given for this example as 0.335 mRad. The NET calculated for this example system is 
0.10 K. Following the reference in Moser this value would be replaced with NET =0.25 K. 
If a FLIR system has been used for a few years, the efficiency of the system will be 
degraded , and the NET will increase slightly. The choice of NET from 0.1 K to 0.25 K is 
reasonable. The resulting computation from the tabulated data for narrow field of view is 
shown in Figure 10. 

B. SELECTED COMPARISON FLIR PARAMETERS FOR WFOV 

Since the range measurements were reported as being made in WFOV, Figure 11 
shows an estimated MRTD for WFOV. In changing the optical field of view, changes will 
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occur in Ax and Ay, the detector angular subtenses, in px through the detector resolution, 
add in the system MTF. To obtain the required WFOV of 15° by 20° a linear magnification 
factor of 3 is required, giving detector angular subtense of Ax and Ay equal to 0.75 mRad. 
Shumaker also provides MDTD in Example 8-16, p. 8-59, with the set of parameters 
which are the same as for MRTD and considered to be appropriate to the Common 

Module FLIR in NFOV mode. The back-end resolution Kb = (rj + ^ + 7 * 5 ^^ in this 

example is also given as 0.335 mRad. NET calculated in this example is also 0.1 ®C. Since 
the range measurements were reported in WFOV, an equivalent MDTD as a function of 
Target Angular Subtense (TAS) for WFOV was also deduced and shown in Figure 12. As 
discussed before, a change of the optical field of view will change Ax and Ay, the detector 
angular subtense, as well as ro, the detector resolution. To obtain the required WFOV, a 
linear magnification factor of 3 is also required, giving detector angular subtense of Ax and 
Ay equal to 0.75 mRad. The detector resolution ro ,which equals Ax, is also 0.75 mRad. 
The back-end resolution in WFOV with Ax=0.75 mRad will be 0.487 mRad. The 
following parameters which are postulated have been discussed above and selected to be 
appropriate to the Common Module FLIR in WFOV. 


MTF(v) calculation for WFOV 


NET 

0.125 “C 

Spatial Frequency ( v ) 

0.75mRad MTF 



(cyc/mRad) 


Ax 

0.75 mRad 

0.00 

1.00 

Ay 

0.75 mRad 

0.25 

0.78 



0.50 

0.55 



0.75 

0.35 



1.00 

0.15 



1.10 

0.085 



1.25 

0.02 



1.30 

0.00 



2.00 

0.00 



2.50 

0.00 



3.00 

0.00 


kr 
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Shumaker [Ref. 8, p. 8-52] shows a recalculation of the system MTF and MRT for a 
decrease of optical magnification, increasing the detector angular subtense to 1 mRad. 
This was taken as a guide in deducing the comparison with Figure 13. An additional check 
point is the expectation of the asymptotic behavior of the MRT(v) at a critical spatial 
frequency between 0.7/Ax and 0.9/Ax. The 1.0 mRad and 0.25 mRad curves in Figure 13 
show the MTF versus spatial frequency for the example system with NFOV and WFOV. 
The cutoff behavior of the MTF gives an indication of the limiting behavior of MRT. The 
estimated MTF (WFOV) is shown in Figure 13 compared with the NFOV and Shumaker’s 
1 mRad model. 

For our comparison MRT, a modification of Shumaker’s MTF computation was 
made with the results shown in Figure 13. The comparison curve lies between the curves 
for the Narrow and Wide Shumaker’s example cases, as expected, with a cutoff very close 
to 1.3 cyc/mRad as appropriate for Ax = 0.75 mRad. 
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VI. DATA COMPUTATIONS AND ANALYSIS 


The process of comparison can be outlined as follows : 

1. The target ship signature Ts was obtained from thermal images obtained with 
the airborne AGA operated by NRaD and compensated to zero range by 
Hughes and McGrath .[Ref 11] 

2. Ship temperature and sea background temperature were converted to in-band 
radiance by application of Planck function. 

3. SEARAD transmission code was used to calculate path transmittance and 
path radiance for the path from the target to the detection /recognition point 
as reported by the operational FLIR operators. 

4. Ship and background radiance were corrected to the detection/recognition 
range using the SEARAD outputs. 

5. Ship and background radiance at detection point were converted to apparent 
temperature using the inverse Planck function. 

6. Apparent temperature difference ATapp at maximum range was compared with 
the postulated MRTD or MDTD respectively calculated for the same target¬ 
sensor locations. If the criteria for detection or recognition are correct, and 
the MRTD and MDTD data are appropriate, the apparent temperature 
difference should match the MRTD or MDTD. 

An example of the computation process is shown in the following section. The 
calculated results of the eight runs for AT and MDTD/MRTD as well as line pairs 
requirement for MRTD criterion at detection range are summarized in Table 10, Table 11, 
Table 12 and Table 13 respectively. The flow charts for the calculation of ATapp, 
MDTD/MRTD and line pairs requirement at detection range are outlined in Figure 6, 
Figure 7, Figure 8 and Figure 9 respectively. 

A. DATA COMPUTATION STEPS 

Run 1, Detection Time : August 04, 04:11 (Zulu), 1992 
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Apparent AT Computation 


Step 1, The data for ship temperature corrected to zero range is taken from 
the Table of Appendix C of Hughes and McGrath [Ref 11] and is reproduced 
here in Table 2 for reference. The average of six ship temperature 
measurements at port aspect is 288.3 “K. The average ship temperature for 
each aspect is shown in Table 9. 


Step 2, The ship temperature is then converted to in-band radiance at zero 
range 


N. 


sgP 

n 


= 


0.95x5.67x10’* x (288.3)% 0-25357^ ^ 300(W / m^Sr) 


where 


e = 0.95, emissivity of R / V POINT SUR. 

= Sjux 288.3 = 023cm xK ^q\ = 1.2227 x 10’’ (Radiation Table) 

X^T = \2^x 288.3 = 0.35cw xK=>q2 = 3.7585 x lO"' (Radiation Table) 

Aq = q2-ql-2.53578x10’ 

The q values have 5 digits as reported in the radiation table (Appendix A) to 
keep the fraction of the radiant existence as close as possible to the integral of 
Planck’s Law. The details of the Proportional Radiation technique are shown 
in Chapter II. [Ref 5] 

Step 3, Run the SEARAD program to obtain transmittance x and path radiance 
Np between FLIR sensor and target. The inputs of the SEARAD code for the 
total 8 runs are tabulated in Appendix E and outputs of the SEARAD 
(transmittances, path radiances and background radiances) are shown in Table 
6 (detection). Table 7 (classification) and Table 8 (identification) separately. In 
Run 1, the transmittance X =0.0369 and path radiance Np=28.7 W/m^Sr. 


Step 4, Ship radiance at full range =A^ t +Ap 

=%'2’+% =30.0x0.0369+28.7 =29.8 (W/m'Sr) 


Step 5, Convert target radiance to equivalent blackbody temperature. 
In-band radiance at full range can be converted into equivalent blackbody 
temperature by use of Planck’s radiation law using an interative computer 



program and look-up table as shown in Appendix B. When Ns+p = 

29.8 WWSr, the = 285.2K. 

• Step 6, From SEARAD output of step 3 , the background radiance Aj at full 

range is 29.6 W/m^ Sr. By conversion with the computer program in Appendix 
B, the corresponding temperature is 284.8 K. 

• Step 7, The apparent temperature difference at full range is then T^+p -Tb =0 4 
^C. 


2. MDTD Computation at Detection Range 

• Step 1, The projected area of the ship on PORT aspect at detection range 
is shown in Figure 4. 

Aj = hv s\n 6+hw cos 6Sin <!)+hi cos0 cos (j) 

= 41.5 •9.75«sin(0.85)+8.8«9.75»cos(0.85)sin(0)+41.5*8.8«cos(0.85)cos(o) 

= 371m^ 

where 

1 = 41.5 m (ship length) 
w = 9.75 m ( ship width) 
h= 8.8 m (ship height) 

<!>= 0.85 “ (elevation angle, sin ' (altitude / R)) 

6=0° (azimuth angle, port aspect) 


• Step 2, The critical dimension of the ship = = 19.3 m 

• Step 3, For pure detection, number of bars N = 1 

• Step 4, The physical dimension to be resolved is x /N=19.3/1=19.3 m 

• Step 5 , Target angular subtense TAS=10^ x A.^ /= 10^ x 371 / (21000p = 
0.84 (mRacP) 

• Step 6 , The typical MDTD curve as a function of Target Angular Subtense 
(TAS) for Shumaker application that is applicable for WFOV is shown in 
Figure 12. From Figure 12 when TAS = 0.8 mRact, the corresponding 
MDTD= 0.2^C 
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3. MRTD Computation 

a. MRTD for Classification 

The observed classification range in Run 1 is 8000 m at elevation angle 0 
==2.2° and azimuth angle (|) =0°. 


• Step 1, For this aspect, the target area is given by 
Ay = Iw^m0+hwcos0s\n(l>+hlcos0cos<f) 

= 41.5 ‘O.TS* sin(2.2)+8.8*9.75* cos(2.2) sin(o)+41.5*8.8‘ cos(2.2) cos(o) 

= 380 

where 

1 = 41.5 m (ship length) 
w = 9.75 m ( ship width) 
h= 8.8 m (ship height) 

9= 2.180 ° (elevation angle, sin ' (altitude / R)) 

^ = 0° (azimuth angle, port aspect) 


• Step 2, The critical dimension of the ship Dc = = 19,5 m 


• Step 3, From the Johnson Criterion, the number of bars N necessary to 
classify the ship is 4 eye x 2 bars/cyc = 8 bars. 

• Step 4, The physical dimension to be resolved is x =D/N=\95I%=2A m 




Step 5, The spatial frequency f characterizing the ship at classification range 
is; 


/.= 


R 

2000x 


8000 

2000*2.4 


(eye / m Rad) 


• Step 6, The MRTD curve as a function of spatial frequency in WFOV for 
Common Module FLIR has been deduced in Chapter V and is shown in Figure 
11. From Figure 11, the asymptotic behavior occurs at spatial frequency equals 
1.33 cyc/mRad. When f=\.l cyc/mRad, the MRTD can not be resolved, 
which indicates that using WFOV is inappropriate in this classification 
application. 
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b. MRTD for Identification : 

The observed identification range in Run 1 is 6000 m at elevation angle 
9 =2.8° and azimuth angle (|) =0° 

• Step 1, The projected area of the ship on PORT aspect for identification is: 

Aj = Iw 9+lm cos 6 ^\n(j)+hl cos 6 cos (j) 

= 41.5 •9.15* sin(2.8)+8.8*9.75‘ cos(2.8) sin(0)+41.5‘8.8* cos(2.8) cos(0) 
= 384 

where 

1 = 41.5 m (ship length) 
w = 9.75 m ( ship width) 
h = 8.8 m (ship height) 

9= 2.8 “ (elevation angle, sin"' (altitude / R)) 

^ = 0° (azimuth angle, port aspect) 

• Step 2, The critical dimension of the ship Dc = = 19.6 m 

• Step 3, From the Johnson Criterion to classify ships, number of bars N 
necessary to classify the ship is 6.4cyc x 2 bar/cyc = 12.8 bars. 

• Step 4, The physical dimension to be resolved is x=T)/A^=19.6/12.8=1.5 m 

• Step 5, The spatial frequency f characterizing the ship at identification range 
is f =6000/(2000x1.5)=2.0 cyc/mRad 

• Step 6, From Figure 11, when f =2.0 cyc/mRad, the corresponding MRTD 
can not be resolved. This may show that using WFOV in this identification 
application is also inappropriate. 

B. LINE PAIRS REQUIREMENT AT DETECTION RANGE 

In Run 1, the calculated apparent temperature difference AT^pp at detection range 
is 0.4 “C. When the target can be resolved at this range, the ATapp must equal MRTD, 
namely MRTD = 0.4 °C. From Figure 11 we know that for WFOV, when MRTD = 0.4 
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°C, the corresponding critical spatial frequency to resolve the target equals 0.9 cyc/mRad. 
From Equation (13) the critical spatial frequency = 0.9 cyc/mRad. At the 

detection range, R = 21000 m, the corresponding number of bars N necessary to resolve 
the target is Dc/%= 1.65 bars, where in Run 1 is 19.3 m. Therefore the line pairs 
require to resolve the target is half of the bars which is 0.83 cycles. 
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Vn. RESULTS AND DATA ANALYSIS 


The calculated apparent temperature difference AT, MDTD and MRTD for the 
experiment are shown in Table 10, Table 11 and Table 12 separately. Line pairs 
requirements to resolve the target at the detection range are shown in Table 13. Data were 
taken during the eight overflights in the period 04:11 to 05:22 (GMT) on August, 1992. 
The first four overflights were flown at approximately 300 m altitude, each pass at a 
different quadrant of the ship. The other four overflights were flown at about 150 m 
altitude. There were three runs approaching the PORT aspect, three runs approaching the 
starboard (STBD) aspect, one run approaching stern (STN) and one run approaching 
BOW aspect. The required meteorological data (including radiosonde measurement) for 
SEARAD code inputs were taken from the meteorological data that were recorded during 
the experiment, and are shown in Table 3 , Table 4 and Table 5. These meteorological 
data were the closest to the aircraft overflights. Ship (RA^ POINT SUR) temperature was 
the average ship temperature taken from the AGA scene corrected for atmospheric effects 
for each aspect (PORT, STBD, STN and BOW). These data were extracted from Hughes 
and McGrath. [Ref. 11] Due to aircraft and balloon flight limitations around the lUV 
POINT SUR, there was at least two hours time difference between the apparent ship 
temperature measurements from the airborne AGA scene and the navy aircraft overflights. 
During the eight overflights, there was no great change in the atmospheric conditions. 

The estimated MRTD curve and MDTD curve for calculating MRTD and MDTD 
at different ranges were derived from Shumaker’s example 8-13 and 8-16 for a generic 
FLIR system in WFOV. The detector angular subtenses Ax and Ay were chosen as 0.75 
mRad to represent a typical FLIR system in WFOV. The NET was chosen as 0.125 K, 
between 0.1 K to 0.25 K, due to degradation of the FLIR sensor after being in service for 
several years. 

In SEARAD, the meteorological data inputs include wind speed (W,), wind 
direction (fVo), visibility (VIS), and air mass character (AM). They also include the sea 
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temperature, aircraft altitude and viewing angle. For analysis purposes, the meteorological 
data were assumed to be the same for each run except for wind speed and wind direction. 
The outputs of SEARAD (atmospheric transmittance t and path radiance Np) at detection 
range, classification range and identification range are shown in Table 10, Table 11 and 
Table 12 separately. Comparisons for each range prediction are shown in the following 
four sections. 

A. DETECTION RANGE COMPARISON 

Table 10 shows Apparent Temperature Difference ATg^pp at the observed 

detection range compared with the MDTD estimated for the postulated system at the same 
range. The data show very large scatter; the variation in the Temperature Difference is 
much more than that of the MDTD for observed ranges from 13 to 22 kilometers. The 
range of discrepancy is from 5 % to 600 % over this small data set. The data used include 
five angles of approach to the target and hence combinations of target temperature and 
apparent size. Comparison between observations with similar aspect, for example runs 5 
with 7 at STBD aspect and 6 with 8 at PORT, show much better internal agreement due 

to approximately the same altitude and viewing angle; However run 2 at 70° azimuth 
shows veiy poor agreement, 

B. CLASSIFICATION RANGE COMPARISON 

Table 11 shows Apparent Temperature Difference and characteristic target Spatial 
Frequency ( fc = Range/(2000 x critical resolution length), see section VI A.3, p.32) at 

the reported classification range. The values shown as MRTD are the values of the 
postulated system WFOV MRTD evaluated at the critical spatial frequency for 
classification, taken as 4 line pairs per critical dimension. The x entries on the MRTD row 
indicate that the WFOV MRTD cannot be evaluated at the required spatial frequency, 
since this falls beyond the asymptotic value. Although the range measurements were all 
reported to be made in WFOV, the comparison suggests that the observed ranges are 
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more appropriate to NFOV operation. The major uncertainty in this table may be the 
FLIR operator's decision that the target was in classification range. In the cases for which 
.comparison is possible, the apparent AT considerably exceeds the postulated wide field 
MRTD. 

C. RECOGNITION RANGE COMPARISON 

Table 12 shows Apparent Temperature Difference and characteristic target spatial 
frequency ( fc = Range / 2000 x critical resolution length ) at the reported 
identification range. The MRTD values are the values of the postulated system in WFOV. 
The MRTD was estimated at the critical spatial frequency for identification, taken as 6.4 
line pairs per critical dimension. The x entries on the MRTD row, as in Table 11, indicate 
that the WFOV MRTD cannot be evaluated at the required spatial frequency. Although 
the range measurements were all reported to be made in WFOV, this comparison still 
suggests that the observed ranges are more appropriate to NFOV operation. The major 
uncertainty in this table may be due to the FLIR operator’s subjective decision that the 
target was in identification range. Again, where comparison is possible, the apparent AT 
considerably exceeds the postulated WFOV MRTD. 

D. RESOLUTION LINE PAIRS AT DETECTION RANGE 

The Johnson Criterion for detection is stated as a required number of resolvable 
line pairs on the critical dimension of the target image. This can be compared with the 
estimated actual line-pairs at the detection range by setting the AT at this range equal to 
MRTD, and evaluating the corresponding spatial frequency. This comparison is shown in 
Table 13 . The required line-pair number is seen to vary from 0.5 to 1.1, compared with the 
Johnson Criterion of 1 The spatial frequency values are read from the MRTD WFOV 
curve for the postulated system, shown in Figure 11. 
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E. DISCUSSION OF THE DATA 


The comparisons of Table 11 through 13 show poor agreement between the 
MDTD/MRTD estimated for the postulated FLIR and the computed target-background 
temperature difference at the appropriate maximum range. No clear systematic 
relationship is apparent between these two quantities, and larger random deviations are 
present. A number of sources can be suggested for the discrepancies. 

The time difference between the apparent ship signature measurement from AGA 
scene with the aircraft overflights measurement may cause some deviation of the 
calculated results. 

The SEARAD code itself, if the inputs were selected properly, should agree to 
within a few degree Celsius with actual sea radiance in the long wave band (8 to 12p.m). 
The results of the SEARAD output (t and Np) and the selection of input parameters may 
also cause the deviation of the computed ATapp. 

The aircrew plays an important role during the measurement. The operator 
generally may not know exactly the target’s characteristics. His primary experience is 
visual ( 0.4 to 0.7 pm ). He can only state if the image on the screen resembles a visual 
scene. In this experiment, the aircrew recorded the detection range when he first saw a 
dot on the FLIR display screen. However, MRTD is a measurement of the performance of 
the entire system which includes the operator, so the measurements of classification and 
identification are subjective, a human interpretation of the displayed information. From the 
data collected by the aircrew in Table 3, the detection ranges for STBD aspect (run 2, 5 
and 7) under similar conditions are 15000 m, 22000 m and 22000 m. The results for run 5 
and run 7 match well with each other, but show great discrepancy with run 2. This 
discrepancy may be due to the different operators’ perception. 

The apparent temperature difference ATapp is a function of the FLIR sensor altitude 
and viewing angle. During some overflights, the sensor altitude and viewing angle changed 
slightly. For example, in run5 and run 7, the sensor altitude was 140 m and 160 m with 
viewing angle 0.36° and 0.43 °, which shows almost the same geometry. The comparison 
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of ATappWith MDTD in Table 10 at STBD aspect show agreement between runs. This is 
also found for run 6 and run 8 at PORT aspect. 

Field of view is also an important consideration for FLIR performance. In normal 
operation a wide field of view ( WFOV ) is initially used to provide situational awareness 
to the operator, but after the target is detected the operator switches to a narrow field of 
view ( NFOV ). During the eight runs in this data set the aircrew reportedly always used 
the WFOV for detection, classification and identification ranges. The results for detection 
shown in Table 10 indicate moderate agreement, but show large discrepancy for MRTD at 
classification and recognition range, as shown in Table 11 and Table 12. The data appear 
more consistent with NFOV measurement. A comparison with the NFOV MRTD might 
show closer agreement. 

Lastly, the “ postulated” FLiR model MDTD and MRTD are based on some 
relatively broad assumptions. The calculated values at detection, classification and 
identification ranges for varies target area and range and meteorological conditions show a 
considerable level of internal consistency, leading to a modest level of confidence in the 
model. 
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Vin. CONCLUSIONS 


The main goals of this thesis were to evaluate the followings : 

• apparent target-background temperature difference (ATapp) at maximum range. 

• MDTD/MRTD of a typical FLIR system applied to PREOS92 experiment. 

• line pairs requirement at detection range to resolve the target. 

In the calculation procedures, the theoretical Planck’s radiation law was applied to 
find ATapp. The propagation code (SEARAD) provided the atmospheric transmittance and 
path radiance. Measured atmospheric parameters were selected as close in time to the 
operational aircraft overflights as possible using PREOS 92 data as inputs to the 
SEARAD code. The ship was required as a rectangular parallelepiped model of the ship’s 
physical length, width and height for MDTD/MRTD evaluation. The geometry data for 
MDTD/MRTD evaluations were also selected from the PREOS 92 experiment 
measurement data. The MDTD/MRTD functions for a generic FLIR in wide field of view 
(WFOV) application were deduced following Shumaker. The Johnson criterion for the 
number of line pairs on a target was employed as a detection criterion and compared to 
the estimated line pairs computed from the data. 

Based on the measured ranges reported in PREOS92 experiment and the calculated 
ATapp at different ranges, the ATapp at the detection range matches better with MDTD 
using the same sensor altitude and viewing angle. However, the measured ranges reported 
are all subjective due to the operator’s interpretation of the displayed information, thus the 
final results of this calculation are all estimations. Furthermore, the calculated MRTD in 
Table 10, Table 11 and Table 12 show that the most of the required spatial frequencies for 
the ranges found in this experiment are higher than the WFOV cutoff spatial frequency 
(1.3 cyc/mRad) , therefore the target should not have been resolved. This suggests that 
using NFOV in this experiment for target classification and identification would be more 
appropriate. 
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Additionally, the calculated line pairs needed at the detection range for MRTD 
criterion show that there is 10% —50 % error range for the 8 runs. Compared to the 
Johnson Criterion quoted as 1 ± 0.25 mRad, this error is acceptable for the calculation 
procedure based on the MRTD curve derived from a theoretical equation used in 
WFOV. This MRTD curve is only a typical curve that is considered to be suitable for this 
application. 

Since only 8 runs that have been applied in this thesis, more runs should improve 
the comparison. 
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Figure 3. Universal Biackbody Curve 

From [Ref. 5] 
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Figure 4. Orientation of Ship Target Model 



N = Number of Bars in the Bar Chart 
R = Slant Range 
TAS = Target Angular Subtense 
= Critical Dimension 
= Physical Dimension to be Resolved 
= Spatial Frequency of Bar 

Figure 5. Relation Between Spatial Frequency, Range R and Physical Dimension % 
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= Ship Temperature at Zero Range 
= Sea Temperature at Zero Range 
r= Transmittance 
Np = Path Radiance 
N, = Background Radiance 
det.r^ = Background Temperature at Detection Range 
= Ship Radiance at Zero Range 
N^^p = Ship Radiance at Detection Range 
= Ship Temperature at Detection Range 

Figure 6. AT Computation Flow Chart 
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Figure 7. MDTD Determination Flow Chart 
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CjTARi^ 



Aj. = Projected Area of the Target 

= Critical Dimension of the Target 

fc = Critical spatial frequency 

L.P= Line pairs required 
to resolve the target 


Figure 8. MRTD Determination Flow Chart 
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Aj. = Projected Area of the Target 

Df. = Critical Dimension of the Target 

= Critical spatial frequency 

N = Number of bars to be resolved 

L.P.= Line pairs required 
to resolve the target 


Figure 9. Line Pairs Calculation Flow Chart 







Spatial Frequency (cyc/mRad) 


Figure 11. MRTD as a Function of Spatial Frequency (WFOV) 




MDTD (C) 



0 0.5 1 1.5 2 2.5 3 3.5 

Target Angular Subtense (mRad)''2 


Figure 12. MDTD as a Function of Target Angular Subtense (WFOV) 





MTF as a Function of Spatial Frequency 



0-5 1 1,5 2 2.5 3 3.5 


Spatial Frequency (Cycles/mRad) 


Figure 13. Comparison of System MTF in WFOV and NFOV 





Table 1. Johnson Criterion 


Discrimination Level 

Cycles Across Minimum Dimension 

Detection 

1.0 

Recognition 

4.0 

Identification 

6.4 


Table 2. Ship Temperature Corrected To Zero Range 
From Thermal Images 

After[Ref. 11] 


Date 

Time (GMT) 

AGA Scene 

Tlshipl'C 





8/04 

01:37:45 

PORT 

14.83 


01:37:47 

PORT 

14.56 


01:39:55 

STBD 

15.06 


01:39:57 

STBD 

15.17 


01:40:01 

STBD 

15.57 


01:46:18 

STBD 

16.68 


01:46:20 

STBD 

17.13 


01:47:49 

PORT 

15.72 


01:47:51 

PORT 

15.97 


01:54:49 

STN 

16.70 


01:56:37 

STBD 

15.72 


01:59:35 

STBD 

15.52 


01:59:37 

STBD 

15.97 


02:00:16 

STN 

16.91 


02:00:51 

PORT 

15.77 


02:01:31 

BOW 

14.79 


02:02:08 

STBD 

16.21 


02:05:50 

BOW 

14.79 


02:06:37 

STBD 

17.19 


02:07:20 

STN 

17.82 


02:07:53 

PORT 

15.18 


02:08:22 

BOW 

14.99 
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Table 3. Data From Aircraft Overflights Measurements 


Aircraft Run 
Numbers 

1 

2 

3 

4 

5 

6 

7 

8 

Run Time (GMT) 

04:11 

04:25 

04:32 

04:41 

04:49 

05:03 

05:11 

05:22 

Ship Aspect 

PORT 

STBD 

STN 

BOW 

STBD 

PORT 

STBD 

PORT 

Aircraft Altitude 
(m) 

313 

307 

310 

309 

139 

155 

164 

162 

Observed 
Detection Range 
(m) 

21,000 

15,000 

13,000 

15,000 

22,000 

15,000 

22,000 

15,000 

Elevation Angle 
for Observed 
Detection Range 
(degree) 

0.8 

1.2 

1.4 

■ 

0.4 

0.6 

0.4 

0.6 

Observed 
Classification 
Range (m) 

8,000 

4,500 

1,800 

8,000 

15,000 

4,600 

9,000 

8,000 

Elevation Angle 
for Observed 
Classification 
Range (degree) 

2.2 

3.8 

9.8 

2.2 

0.5 

1.9 

1.0 

1.1 

Observed 
Recognition 
Range (m) 

6,400 

3,600 

9000 

3,600 

3,600 

2,700 

3,600 

5,500 

Elevation Angle 
for Recognition 
Range (degree) 

2.8 

4.8 

20.0 

4.8 

2.2 

3.2 

2.6 

1.7 
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Table 4. Radiosonde measurement 


4 August 1992 0401 (GMT) 
29 layers_ 


H(kni) 

P(mb) 

T(C) 

RH(%) 

Sets Units for H,P,T,RH 

0.007 

1013.600 

12.200 


ABH 

0.025 

1011.625 

12.280 

HESQH 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 
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Table 5. Meteorological Data And Sea Temperatures For SEARAD Inputs 


Aircraft Run 
Numbers 

1 

2 

3 


5 

6 

7 

8 

Run Time 
(GMT) 

04:11 

04:25 

04:32 

04:41 

04:49 

05:03 

05:11 

05:22 

Air Mass 
Character 

1.6 

1.6 

1.6 

1.6 

1.6 

1.6 

1.6 

1.6 

Visibility (VIS) 

27.08 

27.08 

27.08 

27.08 


27.08 

27.08 

27.08 

Wind Speed Vr 
(m/s) 

10.19 

11.68 

12.12 

13.26 


8.28 

8.05 

7.94 

Wind Direction 
Vd (degree) 

298.911 

289.438 

283.893 

275.374 

277.645 



293.928 

Sea 

Temperature 

T.ea(°C) 

11.7 

11.7 

11.7 

11.7 

11.7 

11.7 

11.7 

11.7 


Table 6. Atmospheric Transmittance (t) And Path Radiance (Np) for Detection 

Range From SEARAD Output File 


Aircraft Run 
Numbers 

1 

2 

3 

4 

5 

6 

7 

8 

Run Time 
(GMT) 

04:11 

04:25 

04:32 

04:41 

04:49 

05:03 

05:11 

05:22 

Transmittance t 

0.0369 

0.0977 

0.130 

0.0965 

0.0092 

0.0672 

0.0150 

0.0689 

Path Radiance Np 
(W/m=Sr) 

28.7 

26.8 

25.8 

26.8 

29.3 

27.7 

29.2 

27.6 

Radiance at 
Observed 
Detection 
Range (w/m^ Sr) 

29.6 

29.2 

29.1 

29.2 


29.3 

29.6 

29.3 

Ship Aspect 


STBD 

STN 

BOW 


PORT 

STBD 

PORT 
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Table 7. Atmospheric Transmittance (t) and Path Radiance (Np) for Classification 

Range from SEARAD Output File 


Aircraft Run 
Numbers 

1 

2 

3 

4 

5 

6 

m 

8 

Run Time (GMT) 

04:11 

04:25 

04:32 

04:41 

04:49 

05:03 

05:11 

05:22 

Transmittance t 

0.261 

0.634 

0.685 

0.252 

0.0726 

0.428 

0.208 

0.241 

Path Radiance Np 
(W/m^Sr) 

21.9 

10.6 

9.0 

22.0 

27.4 

16.8 

23.5 

22.5 

Radiance at 
Observed Classifica 
tion Range 
(W/m^ Sr) 

28..5 

27.3 

27.1 

28.5 

29.2 

27.6 

28.6 

28.4 

Ship Aspect 

PORT 

STBD 

STN 

BOW 

STBD 

PORT 

STBD 

PORT 


Table 8. Atmospheric Transmittance (t) and Path Radiance (Np) for Identification 
Range from SEARAD Output File 


Aircraft Run 
Numbers 

1 

2 

3 

4 

5 

6 

7 

8 

Run Time 
(GMT) 

04:11 

04:25 

04:32 

04:41 

04:49 

05:03 

05:11 

05:22 

Transmittance 

X 

0.341 

0.510 

0.805 

0.507 

0.487 

0.577 

0.498 

0.371 

Path Radiance 
N„ (W/m^Sr) 

19.4 

14.3 

5.5 

14.4 

14.9 

12.3 

14.7 

18.5 

Background 
Radiance at 
Observed 
Identifica tion 
Range 
(W/m'^Sr) 

28.2 

27.5 

28.2 

27.7 

27.4 

26.9 

27.3 

27.8 

Ship Aspect 

PORT 

STBD 

STN 

BOW 

STBD 

PORT 

STBD 

PORT 
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Table 9. Average Ship Temperature (®C) for Each Measurement 

After [Ref. 11] 


Date 

PORT 

STBD 

STN 

BOW 

04 August 92 

153 

16.0 

17.2 

15.1 


Table 10. Comparison of Estimated MDTD with Calculated AT^pp at Detection 

Range 


Run 

1 

2 

3 

4 

5 

6 

7 

8 

Observed 
Range (m) 

21000 

15000 

13000 

15000 

22000 

15000 

22000 

15000 

Aspect 

PORT 

STBD 

STN 

BOW 

STBD 

PORT 

STBD 

PORT 

Af„„CC) 

0.40 

0.98 

1.40 

0.93 

0.14 

0.44 

0.20 

0.40 

MDTDCC) 

0.23 

0.20 

0.20 

0.27 

0.21 

0.14 

0.21 

0.14 


Note : Run 2 is at 70° aspect angle 


Table 11. Comparison of Estimated MRTD with Calculated AT,pp at Classification 

Range 


Run 

1 

2 

3 

4 

5 

6 

7 

8 

Aspect 

PORT 

STBD 

STN 

BOW 

STBD 

PORT 

STBD 

PORT 

Spatial Freq. 
(cyc/mRad) 

m 

1.2 

0.6 

3.3 

3.0 

0.9 

1.9 

1.7 

AT„„ (°C) 

2.1 

4.6 

5.6 

1.8 

0.8 

3.6 

1.9 

2.3 

MRTD("C) 

X 

3.4 

0.2 

X 

X 

0.6 

X 

X 


Note : X indicates MRTD can not be calculated 


Table 12. Comparison of Estimated MRTD with Calculated AT,pp at Identification 

Range 


Run 

1 

2 

3 

4 

5 

6 

7 

8 1 

Aspect 

PORT 

STBD 

STN 

BOW 

STBD 

PORT 

■aiiaf 


Spatial Freq. 
Cyc/mRad) 

2.1 

1.5 

0.4 

2.1 

1.2 

0.9 

IH 


AT,„„ CO 

2.6 

4.1 

3.9 

2.2 

3.9 

4.9 



MRTD (°C) 

X 

X 

0,2 

X 

3.2 

0.6 

3.2 

X 
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Table 13. Line Pairs Requirement at Detection Range For MRTD Criterion 


Run 

1 

2 

3 

4 

5 

6 

7 

8 

Aspect 

PORT 

STBD 

STN 

BOW 

STBD 

PORT 

STBD 

PORT 

MRTD 

at Det. 
Range 

fc (cyc/mRad) 

0.9 

1.0 

1.1 

1.1 

0.3 

0.9 

0.4 

0.9 

Line Pairs 
(cycles) 

0.8 

1.0 

0.8 

0.7 

0.5 

1.1 

0.5 

1.1 


Note : Spatial Frequency is taken from Figure 11 in WFOV 
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APPENDIX A. PROPORTIONAL RADIATION TABLE 


From [Ref. 5] 

The “Proportional Radiation”, q=f(X , T) represents the fraction of the radiant 
exitance which is emitted by a blackbody at temperature T at all wavelengths up to the 
selected values of A,. It is obtained from the integral over the Planck Law for wavelengths 
up to the selected values divided by the integral over all wavelengths, i.e. 


Im(Xt) dX 

< 7 =- 5 ^- 

ffi(XT) dX 


cm. K 

q 

cm . F 

o 

o 

ut 

1.3657,10 ’ 

0-110 

0.0 5 P. 

3■6788.10 ’ 

0 . 150 

0 . 05 ^ 

9 . 17-1 9 . ) 0 ° 

0-160 

0.056 

7 . 1358 . 10 '' 

0 - 170 

0.058 

4.67-15.10® 

0-180 

0.060 

9.6798.10 * 

0.190 

0.06P 

1.9069.10’ 

0.200 

0.06^ 

3.5907.10' 

0.210 

0.066 

6.4902. 10 ' 

0.220 

0-068 

1.1302.10* 

0.230 

0-070 

1.9025.10 * 

0.210 

0 - 07 P 

3.1045.10* 

0.250 

0.071 

4.9236, 10 * 

0.260 

0.076 

7.6070. 10 * 

0.270 

0.078 

1.1473.10* 

0.280 

0.080 

1.6923.10 * 

0 - 290 

0-002 

2.4453. 10 * 

0.300 

0.081 

3 . 4668. 10 * 

0.3 .10 

0.086 

4.8287.10* 

0.320 

0.088 

6.6159.10 * 

.0-330 

0-090 

8.9269.10* 

'0.310 

0.092 

1.187 4.10 '- 

0.350 

0.091 

1.5586. lo '- 

0.360 

0.096 

2.0204 . lo '- 

0.370 

0.098 

2.5885. 10 '• 

0.380 

0-100 

3 . 2.804 . lo '- 

0.390 

0 . 1 1 0 

9.2957, lo '- 

0.100 

0.120 

2.1727.10 ’ 

0.120 

0.130 

4.3066.10’ 

0.110 

Adapted from 

’‘Optoelectronics; 


• cLum • upcoe i ec t ron ics ; 

, McGrAw-Mill Dock Compatiy^ 


X . T 



cm - 

K ^ 

2 . 9053 . 10‘* 

0.160 

5.8057.iO' 

1.3023.10'^ 

0.180 

6.0880-10' 

1.9962. 10"* 

0.500 

6.3191.10' 

2.8858. 10'^ 

0.520 

6.5912.10' 

3.9754. 10'* 

0.510 

6-8146.10' 

5.2613. lO"* 

0.560 

7.0209.10' 

6.733 1. 10'^ 

0.580 

7 . Mi6 . 10‘ 

8.3750. 10'^ 

0.600 

7.3877.10' 

1.0168. lO"' 

0.620 

7.5505.10' 

1.2091. 10"' 

0,660 

7.S102.10' 

1.4122. 10"' 

0.700 

6.0885.10' 

1.6239. 10"' 

0.710 

8.3020.10' 

1.84 23. 10"' 

0.780 

8.1861.10' 

2.0653 . 10"' 

0.820 

8.6155.10‘ 

2.29 11.10"' 

0.860 

8.7810.10’ 

2.5183. 10"' 

0.900 

6 . 9018. 10' 

2.7454. 10"' 

0.910 

9.010H .10' 

2.97 12 . l6"* 

0.980 

9 - 1033. 10' 

3 . 194 7. 10“' 

1 - 00 

9. 1155. 10' 

3.4150. 10"' 

1.10 

9. 3217- 10* 

3.63 14.10"' 

1.20 

9.1532.10' 

3 . 84 3 2 . 10”' 

1.30 

9.5331.10' 

4.0502 . 10"' 

1.10 

9.6301. 10' 

4.2518. 10"' 

1.50 

9.6909 , lO" 

4.44 79.10*' 

1.60 

9.7390. 10' 

4.6382. 10"' 

1.70 

9.7777. 10* 

4.82 27. 10*' 

1.60 

9.809 1 - 10*' 

5.1 7 38.10** 

1-90 

9.6319. 10*' 

5.5012. 10*' 

2,00 

9.8563. 1 0 

7 € o r y and P r a 
•78 . 

ctice " 

A . Cliappe 1 1 . 
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APPENDIX B. LOOK-UP TABLE 


Temperature and Radiance 
Conversion And Look-up Table 

% Numerical Integration for Planck's law 
% Integration by midpoint rule 
% Matlab program with look-up table outputs 

clear; 

Tniin=283.; % <— Guess T to obtain the corresponding radiance 

dT=0.01; 

n=200; 

dispC Temperature Radiance'); 

dispC.'); 

h=6.626e-34; c=2.998e8; k=1.382e-23; 

a=8e-6; b=12e-6; dx=0.001e-6; 

for T=Tmin:dT:Tmin+n*dT 

sf=0.0; 

for x=a:dx:b-dx 

xm=x+dx/2; % midpoint of f(x) 

fl =(2*h* c'^2)/(xm^5); 

12= 1 /(exp((h* c)/(xm*k * T))-1); 

fx=fl*f2; 

sf=sf+fx*dx; 

end; 

disp([T,sf]); 

end; 

end; 
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Look-up Table for Temperature and Radiance Conversions 


Temperature Radiance 


(K) 

(W/m^2Sr) 

283.0000 

28.8410 

283.0500 

28.8673 

283.1000 

28.8937 

283.1500 

28.9201 

283.2000 

28.9464 

283.2500 

28.9728 

283.3000 

28.9993 

283.3500 

29.0257 

283.4000 

29.0521 

283.4500 

29.0786 

283.5000 

29.1051 

283.5500 

29.1316 

283.6000 

29.1581 

283.6500 

29.1846 

283.7000 

29.2111 

283.7500 

29.2377 

283.8000 

29.2642 

283.8500 

29.2908 

283.9000 

29.3174 

283.9500 

29.3440 

284.0000 

29.3706 

284.0500 

29.3973 

284.1000 

29.4239 

284.1500 

29.4506 

284.2000 

29.4773 

284.2500 

29.5040 

284.3000 

29.5307 

284.3500 

29.5574 

284.4000 

29.5842 

284.4500 

29.6110 

284.5000 

29.6377 

284.5500 

29.6645 

284.6000 

29.6913 

284.6500 

29.7182 

284.7000 

29.7450 

284.7500 

29.7719 
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284.8000 29.7987 
284.8500 29.8256 
284.9000 29.8525 
284.9500 29.8794 
285.0000 29.9064 
285.0500 29.9333 
285.1000 29.9603 
285.1500 29.9872 
285.2000 30.0142 
285.2500 30.0412 
285.3000 30.0683 
285.3500 30.0953 
285.4000 30.1224 
285.4500 30.1494 
285.5000 30.1765 
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APPENDIX C. AUXILIARY PROGRAM ( SEARAD CODE ) 


A. INTRODUCTION 

“SEARAD” is a modified M0DTRAN2 code that predicts sea radiance between 
52.63 cm * and 25000 cm'*. MODTRAN is a computer code designed to determine 
atmospheric transmission and radiance at moderate resolution from 0 to 50,000 cm'* . It is 
based on AFGL’s LOWTRAN7 code. SEARAD is based on the Cox-Munk (Cox and 
Munk, 1954,1956) statistical model for wind-driven capillary wave facets. An individual 
facet is chosen and a specific slope is given with respect to the local horizon. The facet is 
allowed to reflect the sky and sun and emit thermal black body radiation toward an 
observer. We can find the total radiance by applying proper weight to the facet and 
integrating over all facets within the observer's field of view. 

“SEARAD” can be used in a spectral range from the visible to the far infrared 
regime. It is a DOS-compatible program that runs on a personal computer and calculates 
sea radiance. The operation of SEARAD is very similar to the original MODTRAN2 code 
except for some additional new routines to compute the spectral radiance from sun , sky 
and sea. A new logical switch “SeaSwitch ” is required in the input file to activate the 
radiance, calculation. Sun glint is also include in the sea radiance prediction. Preliminary 
comparisons with data show that SEARAD has an accuracy approximately “ several °C in 
the infrared”. [ Ref 1, p. 4 ] 

SEARAD is currently designed for a single pixel and takes about 10 s to execute. 
Each time a new geometry is chosen by the user. SEARAD recalculates the source 
radiance and the path radiance and transmission. However, only the path properties 
change significantly from one pixel to the next in an ocean image. 

B. UNZIP SEARAD (INSTALLATION) 

The following statement tells how to unzip and run "SEARAD". The zipped files 
span three disks altogether. The last disk will be disk 3, the first disk will be disk 1. These 
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commands assume that the hard disk in the computer is drive c:, and that the 3 1/2 inch 
disks will be loaded into drive a:. If another drive is used instead for the 3 1/2 disks, the 
drive letter (b: for example) should be substituted for a: in the following instructions. 

Sentences in parentheses are physical acts; sentences without parentheses are 
commands that should be typed on the keyboard. 

To unzip: 

c:\ mkdir searad 

cd searad (Insert diskl in drive a; ) 
copy a:\pkunzip.exe 
pkunzip a:\source 
pkunzip a:\input 

pkunzip a:\exe (Follow instructions for disk insertion.) 

The files should all be unzipped at this point. The directory c:\searad should now contain: 
zip utility: pkunzip.exe 

source code : Modnn [nn=10 to 22 ( except mod21.for) 

Input file: TapeSxxx.std 
Log of mods: Note.txt 

Executables : Searad.exe, Dirac, and DOSxmsf exe 
To run Searad : 

copy tapeSrad.std tapeS 

searad 

type out 

This should show the output file corresponding to the input file " TapeSrad.std ". 

C. SEARAD SAMPLE DESCRIPTION 

In this section an example is provided to show how to use SEARAD to predict 
ocean radiance. An input file called " TapeSrad.std " (Table 14) employs a 1976 U.S. 
standard atmosphere to calculate ocean radiance observed at a zenith angle of 100 degree 
fi-om a height of 23 m. The Navy aerosol model is used. The calculation is done with 
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multiple scattering at low spectral resolution (LOWTRAN 7) for a single wave number 
(945 cm) in the long wave band. The meanings of each " Card " field in the input file 
(TapeSrad.std) in SEARAD are interpreted in details in Appendix D. [Ref 1, pp. 3-12 ] 
The following DOS commands will calculate ocean radiance and print results 
copy tapes. std tapeS 
searad 
type out 

These commands produce an output file called “out” (Table 15) . The four 
contributions to ocean radiance (path to footprint, sea emission, sky reflection, and sun 
glint) are listed at the end of the “out” file. “TBOUND” in the input file here is 
interpreted as sea temperature. 

There are two new parameters in the input file at the end of the third line: "90.000" 
and “T”. The “T” (stands for true), which may appear anywhere in columns 76 through 
80 of the third line of the input file (at the end of card 3), is a new logical parameter 
Seaswitch.” The parameter is required or a fatal error will be generated in the input file. 
“Seaswitch” controls the sea radiance calculation. When “Seaswitch” is equal to “T,” the 
sea radiance calculation will be allowed provided certain other conditions are met. When 
it is equal to “F” (stands for false), the sea radiance calculation will be prevented under all 
conditions and the program will execute as originally released by the Air Force. 

The “90.000,” which may appear anywhere in column 66 to 75 of the third line of 
the input file (near the end of the input file), is a new floating parameter. The program will 
run whether the parameter is included in the input file or not. It is the azimuth of the 
upwind direction measured from the line of sight in degrees positive East of North. If the 
field is blank, meaning that the observer is looking directly into the wind. In the input file 
in Table 14, the parameter "90" means that the wind is blowing fi'om right to left, 
perpendicular to the direction of observation. [Ref 1, pp. 4-5] 
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SEARAD MODEL 


D. 

The SEARAD model computes four contributions to sea radiance. The assumption 
for the model is that the strength of interaction between an optical ray and a capillary 
wave facet is given by the facet area projected normal to the ray. It does not include 
multiple reflections, shadowing and gravity waves. It also ignores polarization. Each of the 
four contributions is shown in Figure 14. 

The first contribution is path radiance (Apath), which is shown on top of Figure 14. 
The footprint of a single pixel in an image of the sea is indicated by the wavy line. The 
footprint is observed by a receiver at the end of a ray whose zenith angle at the footprint is 
Or . Apati, represents the spectral radiance in WW /sr/cm along the path fi'om the footprint 
to the receiver. 

The second contribution is reflected sky radiance (Asky). The spectral radiance Ns 
from part of the sky reaches the footprint along a ray whose zenith angle is 6>s. The 
contributions are the summation from all portions of the sky after specular reflection by 
the appropriate facets within the footprint. Asky is the sum leaving the footprint at zenith 
angle ft . During its path to the receiver, Asky is attenuated by the path transmission Tpau,. 

The third contribution is reflected solar radiance ( Asun ), sun glint. The radiance 
No from the solar center arrives at the footprint along a path whose zenith angle is ft,. 
Within the footprint most facets deflect the solar ray away from the receiver and are 
rejected, but some facets are retained because they deflect the ray specularly toward the 
receiver along a path with zenith angle ft. Asun is the spectral radiance leaving the footprint 
after summation over rays arriving from all portions of the solar disk. Nsun is also 
attenuated by the path transmission Tpath ■ 

The fourth contribution is thermal black body emission ( Asea )• In this portion each 
facet emits a spectral radiance Abb given by Planck's equation for a black body. The sea 
temperature is equal to the value of TBOUND in the input file. The spectral emissivity of 
a given facet in the direction of the receiver is specified by the slope of that facet and the 
value of ft. Nsea is the thermal spectral radiance leaving the footprint for the receiver after 
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summation over all facets within the footprint. A^ea is also attenuated by path transmission 
^patii after leaving the footprint. 

The symbol p in Figure 14 represents the reflectivity of sea water. It is calculated 
from Fresnel's equations. The total spectral radiance l^(v) received at wave number 
v( cm-1 ) is given in Equation (19) 

Mv) 

— ^ path 

where f(v) stands for the spectral responsivity of the receiver. 

The design of SEARAD is such that the path (Npath, Tpath) and source (TV,, No, TVbbj 
values are taken from the original MODTRAN 2 while Fresnel reflection (p) and slope 
integrated values (TVsky ,Nsm, NseJ are introduced in new subroutines. Integration of 
Equation(19) over the wave number band specified in the input file ( TapeS ) is 
implemented in a modification of subroutine "TRAN" to produce the band integrated 
values for sea radiance given in the output file (out). [ Ref 1, pp. 6-7 ] 
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Table 14. SEARAD Input File Sample (TapeSrad.std) 

F63 1 1000000000 288.15 0.0 
3 0 0 5 0 0 10.00 10.00 10.00 0.00 0.00 
0.023 0.00 100.00 0.0 0.0 0.0 0 90.0 T 

940 950 10 5 

0 


"TAPE5RAD.STD" for Sky3.in 

Lowtran radiance 

1976 U. S. Standard Atm 

Sea switched on 

Slant Path to space; 

23 m 

100 degrees 

Wind az 90.00 wrt Look 
Thermal mode 
TBOUND = 15 C = Tsea 
No MS reflections: 

SALE = 0. 

Air Mass 5 
Visibility 10 km 
Current wind 10 m s-1 
Average wind 10 m s-1 
Narrow band at 945 cm-1 



Table 15. SEARAD Output File Sample (out) 

***** SEARAD, A MODIFICATION OF LOWTRAN7 ***** 


DATE: 06/09/1996 
THERMAL RADIANCE MODE 
MULTIPLE SCATTERING USED 
MARINE AEROSOL MODEL USED 
WIND SPEED = 10.00 M/SEC 

WIND SPEED = 10.00 M/SEC, 24 HR AVERAGE 

RELATIVE HUMIDITY = 50.00 PERCENT 
AIRMASS CHARACTER = 5.0 

VISIBILITY =10.00 KM 
SLANT PATH TO SPACE 
HI = .023 KM 

HMIN = .000 KM 

ANGLE = 100.000 DEG 

FREQUENCY RANGE 

IVl = 940CM-1 ( 10.64 MICROMETERS) 

IV2 = 950CM-1 ( 10.53 MICROMETERS) 

IDV = 10 CM-1 

IFWHM = 5 CM-1 

IFILTER = 0 

SEA AT 288.15 K REPLACES BLACK BODY BOUNDARY 
UPWIND = 90.000 DEG EAST OF LINE OF SIGHT 

ZERO RANGE VALUES 

SEA EMISSION = .72334 W M-2 SR-1 (AV. EMISS. .7989) 
SKY REFLECTION = .06266 W M-2 SR-1 
SUN GLINT = .00000 W M-2 SR-1 

TOTAL RADIANCE = .78599 W M-2 SR-1 
BLACK BODY TEMP.= 6.7 C 
FULL RANGE VALUES 

PATH TO FOOTPRINT = .01814 W M-2 SR-1(AV. TRANS. .9776) 
SEA EMISSION = .70712 W M-2 SR-1 

SKY REFLECTION = .06125 W M-2 SR-1 
SUN GLINT = .00000 W M-2 SR-1 

TOTAL RADIANCE = .78652 W M-2 SR-1 
BLACK BODY TEMP = 6.7 C 
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APPENDIX D. DESCRIPTION OF EACH CARD IN SEARAD INPUT FILE 


In this input file, the format is very important. Each number must be in the correct 
location. Below is a listing of the input file parameters in SEARAD and the values used in 
evaluating the sea radiance. The definition of each value is included. The xx's serve as 
place holders for actual values. 

CARD 1 (line one ) 

F => use LOWTRAN 7 
7 => MODEL : radiosonde data used 

3 ITYPE ; slant path to space 

1 => lEMSCT : program execution in radiance mode 

1 => EMULT : program executed with multiple scattering 

"M1,M2,M3" are used to modify or supplement the altitude profiles of 
temperature and pressure, water vapor, and ozone "M4,M5,M6" 

XX => Ml : default temperature and pressure to specified MODEL 
atmosphere (Ml=l~6) 

XX => M2 ; default HjO to specified MODEL atmosphere ( M2=l~6 ) 

XX => M3 : default ozone to specified MODEL atmosphere ( M3=l~6 ) 

XX => M4 : default CH4 to specified MODEL atmosphere ( M4=l~6) 

XX M5 ; default N 2 O to specified MODEL atmosphere ( M5=l~6 ) 

XX => M6 : default CO to specified MODEL atmosphere ( M6=l~6 ) 

1 MDEF use default profile for C02, 02, NO, S02, N02, NH3, HN03 

1 => IM ; radiosonde data to be read initially 

1 => NORPT : minimize printing of transmittance 

XX => TBOUND : boundary temperature (®k ) 


75 






0.00 :=> SALE : surface albedo ( 0.0 to 1.0 ) 

CARD 2 (line 2) 

3 => IHAZE ; navy maritime extinction, sets own VIS. 

0 => ISEASN : Summer for model 7 

0 => IVULCN ; default to stratospheric background 

16 => ICSTL ; air mass character for 1.6. (1 for open ocean - 10 for strong continental 

influence) ICSTL now allow non-integer values. Integers from 
1 to 10 inclusive will be taken as is, while integers from 11 to 300 
inclusive will be divided by 10. Integers outside the range of 1 to 
300 will be reset to a default value of "3.0". 

0 => ICLD : no clouds or rain 

0 r:> IVSA : not used 

XX => VIS : meteorological range 

XX => WSS : specifies current wind speed ( m/s) 

XX => WSS : 24 hour average wind speed ( m/s ) 

0.00 => RAINT : rain rate ( mm/hr ) 

0.00 GNDALT altitude of surface relative to sea level (km ) 

CARD 2c (line 3 ) 

33 => ML ; number of atmospheric layers to be inserted 

0 => IRDl : no read 

0 ^ IRD2 ; no read 

text TITLE ; user defined text 
CARD 2c 1 ( next 33 line ) 

XX => ZMDL : altitude of layer ( km ) 

XX P : pressure at layer (mb) 
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, XX T ambient temperature (°C) 

XX nr> JCHAJR.(RH): relative humidity (% ) 

ABH set units for P,T,RH 
CARD 3 (line 37 ) 

XX => HI : initial altitude (FUR sensor altitude) (km) 

0.00 => H2 ; final altitude (seasurface altitude) (km) 

XX => ANGLE : zenith angle (degree) 

0.000 => RANGE ; path length ( km ) 

0.000 => BETA : earth center angle subtented by H! and H2 (degree) 

0.000 => RO ; default radius of earth (,km ) 

0.000 => LEN normal operation of program ( 1 for long path through tangent height) 
90.00 => Psi : azimuth of the upwind direction from line of sight positive East of North 
(degree) 90 “ means that the wind is blowing from right to left, 
perpendicular to the direction of observation 

T ^ True . sea radiance will be carried out ( F means that sea radiance will be 
prevented under all conditions) 

CARD 4 (line 38 ) 

830 => VI : initial wave number ( cm ' ) 

1250 V2 : final wave number ( cm ’) 

10 => DV : wave number increment ( cm ' ) 

0 no filter 

CARDS 

0 = 5 > IRPT : end program ( 3 means read CARD 3 again) 
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APPENDIX E. SEARAD INPUT DATA 


Detection Range, RUN 1 

F7311666 6 6 6 1 11 284.70 

3 0 0 16 0 0 27.080 10.19 10.50 0.000 0.000 


29 0 

0.007 

0 Aug.04 1992,0411 PORT Aspect 
1013.600 12.200 94.0 ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 


0.3130 0.000 90.850 0.000 0.000 0.000 1 298.911 T 

830 1250 10 10 0 

0 


0.000 
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6 


6 


6 


Detection Range, RUN 2 
F 7 3 1 1 6 6 6 

3 0 0 16 0 0 27.080 11.68 10.50 


29 0 

0.007 

0 Aug.04 1992,0425, 
1013.600 12.200 

STBD 

94.0 

ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.3072 

.000 91.197 .000 

.000 1 

830 

1250 10 

10 

0 



0 


1 1 

.000 .000 


1 289.438 T 


284.70 


0.000 
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Detection Range, RUN 3 

7 3 1 1 6 6 6 6 6 6 1 1 1 284.70 0.000 

3 0 0 16 0 0 27.080 12.12 10.50 .000 .000 


29 0 

0.007 

0 Aug.04 1992,0432, 
1013.600 12.200 

STN 

94.0 

ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 


0.3109 .000 91.390 .000 .000 0.000 1 283.893 T 

830 1250 10 10 0 

0 


81 




6 


1 1 284.70 0.000 


Detection Range, RUN 4 


7 3 

1 1 6 

6 6 

6 

6 

3 0 

0 16 0 

0 27.080 13.26 10.50 

29 0 

0.007 

0 Aug. 04 1992, 0441, BOW Aspect 
1013.600 12.200 94.0 ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.3088 

.000 91.209 .000 

.000 0.000 

830 

1250 10 

10 

0 



0 


1 

.000 .000 


1 275.374 T 


82 



Detection Range, RUN 5 


7 3 

1 1 6 

6 6 

6 

6 

3 0 

0 16 0 

0 27.080 12.04 10 

29 0 

0 Aug. 04 1992, 0441 

0449, STBD 

0.007 

1013.600 

12.200 

94.0 

ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12,390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21,700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3,900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10,0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.1393 

.000 90.364 .000 

.000 0.( 

830 

1250 10 

10 

0 



0 


6 1 
.000 .000 


1 277.645 


1 


T 


284.70 




0.000 


83 


Detection Range, RUN 6 

F7311666 6 6 6 1 11 284.70 

3 0 0 16 0 0 27.080 8.277 10.50 .000 .000 

29 0 0 Aug. 04 1992, 0503, PORT Aspect 


0.007 

1013.600 

12.200 

94.0 

ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 


0.1548 .000 90.570 .000 .000 0.000 1 295.427 T 

830 1250 10 10 0 

0 


0.000 


84 




Detection Range, RUN 7 

F7311666 6 6 6 1 11 284,70 

3 0 0 16 0 0 27.080 8.277 10.50 .000 .000 


29 0 
0.007 

0 Aug. 04 1992,0511, STBD Aspect 
1013.600 12.200 94.0 ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.1640 

.000 90.428 .000 

.000 0.000 

830 

1250 10 

10 

0 



0 


0.000 


85 


Detection Range, RUN 8 

F 7311666 6 6 6 1 11 284.70 0.000 

3 0 0 16 0 0 27.080 7.943 10.50 .000 .000 


29 0 
0.007 

0 Aug. 04 1992, 0522, PORT 
1013.600 12.200 94.0 ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 


0.1615 .000 90.595 .000 .000 0.000 1 293.928 T 

830 1250 10 10 0 

0 


86 


1 1 284.70 0.000 


Classification, RUN 1 


7 3 

1 1 6 

6 6 

6 

6 

3 0 

0 16 0 

0 27.080 10.19 10.50 

29 0 

0.025 

0 Aug. 04 1992, 0411 PORT Aspect 
1011.625 12.280 92.4 ABH 

0.007 

1013.600 

12.200 

94.0 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10,400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.3130 

.000 92.180 .000 

.000 0.000 

830 

1250 10 

10 

0 



0 


6 1 
.000 .000 


1 298.991 T 


87 


6 


6 


6 


1 


I 1 284.70 0.000 


Classification Range, RUN 2 
>,F .7 3 1 1 6 6 6 

3 0 0 16 0 0 27.080 11.68 10.50 

29 0 0 Aug.04 1992,0425, STBD Aspect 


0.007 

1013.600 

12.200 

94.0 

ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.3072 

.000 93.853 .000 

.000 1 

830 

1250 10 

10 

0 



0 


.000 .000 


1 289.438 T 


88 



6 


1 


1 1 284.70 0.000 


Classification Range, RUN 3 


7 3 

1 1 6 

6 6 

6 

6 

3 0 

0 16 0 

0 27.080 12.12 10. 

29 0 

0.007 

0 Aug.04 1992,0432, STN Aspect 
1013.600 12.200 94.0 ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88,0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.3109 

.000 97.788 .000 

.000 0.0( 

830 

1250 10 

10 

0 



0 


.000 .000 


1 283.893 T 


89 




Classification Range, RUN 4 


7 3 

1 1 6 

6 6 

6 

6 

3 0 

0 16 0 

0 27.080 13.26 10 

29 0 

0 Aug. 04 1992, 0441, BOW 

Aspect 

0.007 

1013.600 

12.200 

94.0 

ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.3088 

.000 92.150 .000 

1 .000 0.1 

830 

1250 10 

10 

0 



6 

.000 


1 

.000 


1 1 284.70 


1 275.374 


0.000 


90 



1 I 284.70 0.000 


Classification Range, RUN 5 


7 3 

1 1 6 

6 6 

6 

6 

3 0 

0 16 0 

0 27.080 12.037 10.5i 

29 0 

0 Aug. 04 1992, 0441, 0449, STBD 
Aspect 

0.007 

1013.600 

12.200 

94.0 

ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.1393 

.000 90.545 .000 

.000 0.000 

830 

1250 10 

10 

0 



0 


6 1 
0.000 .0.000 


1 289.911 T 


91 



Classification Range, RUN 6 

F7311666 6 6 6 1 11 284.70 0.000 

3 0 0 16 0 0 27.080 8.277 10.50 .000 .000 

29 0 0 Aug. 04 1992,0503, PORT Aspect 


0.007 

1013.600 

12.200 

94.0 

ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.1548 

.000 

91.940 .000 .000 

830 

1250 

10 10 

0 



0 


92 



6 


6 


6 


1 


1 1 284.70 0.000 


Classification Range, RUN 7 
<F 7 3 1 1 6 6 6 

3 0 0 16 0 0 27.080 8.049 10.50 


29 0 

0.007 

0 Aug. 04 
1013.600 

1992,0511, STBD 
12.200 94.0 

Aspect 

ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 


0.1640 .000 91.027 .000 .000 0.000 

830 1250 10 10 0 

0 


.000 .000 


1 289.911 T 


93 



6 


6 


6 


1 


1 1 284.70 0.000 


Classification Range, RUN 8 
F 7 3 1 1 6 6 6 

3 0 0 16 0 0 27.080 7.943 10.50 


29 0 

0.007 

0 Aug. 04 1992, 0522, PORT Aspect 
1013.600 12.200 94.0 ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.1615 

.000 91.125 .000 

.000 0.000 

830 

1250 10 

10 

0 



0 


.000 .000 


1 293.928 T 


94 




Identification Range, RUN 1 

,F7311666 6 6 6 1 11 284.70 

3 0 0 16 0 0 27.080 10.19 10.50 0.000 0.000 


29 0 

0.007 

0 Aug.04 1992,0411 PORT Aspect 
1013.600 12.200 94.0 ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

: 0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 


1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5,0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19,100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 


0.3130 0.000 92.803 0.000 0.000 0.000 1 298.991 T 

830 1250 10 10 0 

0 


0.000 


95 




Identification Range, RUN 2 

F7311666 6 6 6 1 

3 0 0 16 0 0 27.080 11.68 10.50 .000 .000 

29 0 0 Aug.04 1992,0425, STBD Aspect 


1 1 284.70 


0.007 

1013.600 

12.200 

94.0 

ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.3072 

.000 94.818 .000 

1 .000 1 

830 

1250 10 

10 

0 



0.000 1 289.438 T 


0.000 


96 



6 


1 


1 


1 


284.70 


Identification Range, RUN 3 


7 3 

1 1 6 

6 6 

6 

6 

3 0 

0 16 0 

0 27.080 12.12 10.50 

29 0 

0.007 

0 Aug.04 1992,0432, STN Aspect 
1013.600 12.200 94.0 ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.3109 

.000 109.880 .000 

.000 0.000 

830 

1250 10 

10 

0 



0 




.000 .000 


1 283.893 T 


97 



98 




Identification Range, RUN 5 

F731I666 6 6 6 1 

3 0 0 16 0 0 27.080 12.04 10.50 .000 .000 


29 0 
0.007 

0 Aug. 04 1992, 0441, 0449, STBD 
1013.600 12.200 94.0 ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.1393 

.000 92.182 .000 

.000 0.000 

830 

1250 10 

10 

0 



284.70 0.000 


99 



6 


1 


1 1 284.70 0.000 


Identification Range, RUN 6 

•F .7 3 1 1 6 6 6 6 6 

3 0 0 16 0 0 27.080 8.277 10.50 .000 

29 0 0 Aug. 04 1992, 0503, PORT Aspect 


0.007 

1013.600 

12.200 

94.0 

ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 

6.0 

ABH 

0.771 

926.800 

25.600 

4.0 

ABH 

1.525 

850.000 

21.700 

7.0 

ABH 

3.174 

700.000 

12.400 

5.0 

ABH 

3.182 

699.300 

12.300 

5.0 

ABH 

4.017 

632.100 

5.200 

14.0 

ABH 

4.514 

594.600 

3.900 

4.0 

ABH 

4.974 

561.600 

0.000 

6.0 

ABH 

4.984 

560.900 

0.000 

5.0 

ABH 

5.891 

500.000 

-6.700 

7.0 

ABH 

6.295 

474.700 

-10.400 

11.0 

ABH 

6.856 

441.000 

-13.700 

6.0 

ABH 

6.880 

439.600 

-13.800 

25.0 

ABH 

6.912 

437.800 

-14.100 

8.0 

ABH 

7.590 

400.000 

-19.100 

4.0 

ABH 

8.037 

376.600 

-21.800 

5.0 

ABH 

8.184 

369.100 

-23.100 

10.0 

ABH 

8.232 

366.700 

-23.500 

36.0 

ABH 

8.284 

364.100 

-24.000 

14.0 

ABH 

0.1548 

.000 

93.235 .000 .000 

830 

1250 

10 10 

0 



0 


.000 


295.427 T 


100 



Identification Range, RUN 7 

^^7311666 6 6 6 1 11 284.70 

3 0 0 16 0 0 27.080 8.049 10.50 .000 .000 


29 0 
0.007 

0 Aug. 04 1992,0511, STBD Aspect 
1013.600 12.200 94.0 ABH 

0.025 

1011.625 

12.280 

92.4 

ABH 

0.050 

1008.881 

12.390 

90.2 

ABH 

0.075 

1005.709 

12.499 

88.0 

ABH 

0.120 

1000.000 

12.700 

84.0 

ABH 

0.208 

989.500 

12.300 

84.0 

ABH 

0.253 

984.200 

12.200 

84.0 

ABH 

0.270 

982.300 

12.900 

83.0 

ABH 

0.365 

971.300 

16.700 

30.0 

ABH 

0.543 

951.400 

23.700 

6.0 

ABH 

0.558 

949.800 

24.300 
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